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Charge asymmetries for the charm meson production (D* + -D*~ , D*°~D*° and 
D+-D~) have been estimated for the COMPASS kinematic conditions in the frame- 
. work of perturbative recombination model. Mass corrections have been taken into 

o 

© : account in the calculations. The large asymmetry for Df^Dj production has been 

£^ ' predicted. 

m 

1. INTRODUCTION 

> 

A production of hadrons is a good testbed for the research of quark and gluon interactions at 
^ ■ the distance of color confinement. 



in 
o 



On the one hand the intermediate value of c-quark mass complicates the process description in 



' the QCD frame work, on the other hand this quark mass value allows us to research heavy quark 
interaction at the distances, which are larger than there Compton wave length. 

The production asymmetry of charmed particles and antiparticles has been observed in many 
high-energy experiment. For the first time such asymmetry had been discovered in the charm 
hadronic production [1]. The asymmetry value in that process depends on the quark structure of 
initial hadrons, and can be naturally explained by taking into account the interaction with the 
hadronic remnant [2]. 

It is more amazing that a such charge asymmetry exists in the charm photoproduction in the 
photon fragmentation region [3-5]. E691 [3] and E687 [5] Collaborations observed the statistically 
significant asymmetry value for D + — D~ , D° — D° and D* + — D*~ yields. This fact speaks well for 
the essential role of valence quarks of the initial hadron even in the region of photon fragmentation 
at their photon energy. 

The observed asymmetry cannot be explained in the frame work of pQCD and factorization 



* Electronic address: aber@ttk.ru 
^Electronic address: likhoded@mx.ihep.su 



2 



theorem for the inclusive D-meson spectra. It is well known, that in the leading order on a s of 
perturbative theory D- and D-meson one gets identical spectra. At NLO calculations [6] it is 
appeared a small asymmetry due to the quark-photon interaction: 



Indeed, the NLO-calculations for the photon energy of E 1 = 200 GeV leads to the asymmetry 
value which is smaller than the experimental one by the order of magnitude. Thus, the hard pro- 
duction of cc-pair followed by the fragmentation into hadrons can not describe a such asymmetry. 
It is worth to note that this situation is analogous to the situation with weak decays of D-meson, 
where spectator mechanism can not explain lifetime difference between charged and neutral D- 
meson. So, one need to take into account the interaction with a charmless component of meson, 
as well the interference terms. 

Besides of an experimentally observed asymmetry there is another problem in using the fac- 
torization approach to estimate D-meson production cross section. The thing is that at small 
transverse momenta one can not represent D-meson production as cc-pair production with further 
c-quark fragmentation into D-meson, because there is an interference between the light quark pro- 
duced during c-quark fragmentation and that from a hadronic remnant. That taking into account 
the interference contribution one could determine the region, where the factorization approach is 
valid. 

We have tried to explain in [7] the both discussed effects (the asymmetry and the interference) 
by including into perturbative calculation in addition to heavy quarks the light ones needed to 
obtain a hadron with desired quantum numbers. In so doing the light quark has a mass m q ~ 
Aqqq ~ 300 MeV, which plays a role of an infrared cut-off. However, the inclusion of the 
light one essentially complicates calculations. For example, the production process in gluon-gluon 
interaction is described by 36 diagrams of 0(a A s ) order, the production process in photon-photon 
interaction is described by 24 diagrams of 0(aa^) order. 

The principal features of the discussed model are as follows: 

1. At large transverse momenta the model predictions coincide with the predictions of the 
fragmentation mechanism. The cross section dependence on D-meson transverse momentum 
P? have the standard factorized form: 



7 + g 



c + c + q. 



(1) 



da(ij -> D + X) 



J d>P± 



f da(ij — > cc) 



fc^D(z D )S [z D r \ dz D dp c ± 




(2) 
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were i and j are interacting partons, p\ is a transverse momentum of c-quark, and da(ij — > 
cc)/dpc is the production cross section of cc-pair in the leading order. 

2. In small transverse momenta region there is appeared a new contribution that leads to power 

2 

TfL 

corrections — This contribution violates the factorized equation (2). 
Vt 

3. At small pt and large longitudinal momenta (i.e. at large x) it is appeared a strong inter- 
ference and the production cross section can not be described by (2). 

Note that the fragmentation function / c -kd(z.d) in (2) is calculated perturbatively in the frame 
work of the discussed model and has the form: 



ag(Q( 3 Si)) , 
Jc^dM = 5 1 {z D ,r), 



m 2 n m c 

q 



(3) 



where r = m q /m c and / is a function on Zr> and r. Let us to note, that the evaluated function is 
close in form to Peterson fragmentation function [10]. 1 

Variables (0( 3 Si)) and ((O 1 ^)) in the relativistic limit correspond to a wave function squared 
at the coordinate origin. In the discussed model these variables are considered as free parameters. 
The values of {0( 3 Si)) and {(O 1 S )) have been chosen to satisfy the normalization condition 



♦ DM = J fc->DM dz D , (4) 



where W C _^ D (*) is the probability of fragmentation c — > taken from the experimental date on 
e + e~-annihilation. As it is seen from (3) and (4), the values of m q , a s and (0( 3 Si)) (or (O( 1 S , ))) 
are correlated to each other. For example for m q = 0.3 GeV and a s = 0.3 GeV one gets 

(0( 3 S 1 )) ~ (OCSJ) = 0.25 GeV 3 . 

These values allow to describe the experimental data on photo- and electroproduction at HERA 
collider [8]. 

It is worth to emphasize the role of the following process for the charm production: 

7 + q — D^\cq) + c. (5) 



1 It is important to stress, that the common belief in nonperturbative nature of Peterson fragmentation function 
is not valid. The dependence on z in that function is caused by the heavy quark propagator. It had been shown 
in the original work [10], that to obtain this dependence one need to take into account only the right pole of 
propagator l/(E — E') in the infinite reference frame. 
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Due to q <-> q asymmetry in the structure functions of the initial hadrons the process (5) leads to 
a charge asymmetry in the -D^-meson yields. This process dominates at low energies and large 
values of Feynman variable Xp — -7=- where is a longitudinal momentum of D-meson. At high 
energies at small xf the contribution of such process is negligible. At kinematic conditions of 
HERA collider the process (5) contributes only several percents into the total cross section value, 
which is within experimental errors [7]. 

Following [9] we refer the process (5) as a perturbative recombination. The detailed analysis 
of the perturbative recombination reveals that the cross section behavior depends rather weakly 
on a light quark mass. The analytic calculations in the limit as x q — > 0, where x q is a momentum 
fraction of D-meson carried by the light quark, exhibit several interesting features of the process 
(5). That, accordingly to the results of [9], in the initial light quark direction {6 = 0) 



da (7 + q — > (cq) + c) 



256tt 

a s (6) 



da{j + q — > c + c) 

for both (cq)( l s ) and (cg)( 3 Si)-states. One can see that in spite of an additional power of a s the 
cross section value of the process (5) is comparable to the cross section value of the photon-gluon 
fusion. In the photon direction (6 = ir) the production cross section is suppressed by additional 
factors of m 2 c / s [9]: 

^a s ^ for % 

(7) 



da (7 + q — > (cq) + c) 



da(7 + g — > c + c) 



fa, 1 ? for 3 ^ 



At large transverse momenta p± the cross section is suppressed by an additional factor of m 2 c /p\ 

and falls of like . 

PI 

2. CALCULATION RESULTS 

The cross section dependence on scattering angle for perturbative recombination is shown in 
Fig. lfor both vector and pseudoscalar states of (cg)-system. For the sake of comparison the 
cross section behaviour is presented for the case of limit x q — > [9], as well as for the nonzero 
mass of light quark m q = 0.3 GeV. One can see that at low energies the production cross section 
for pseudoscalar (cg)-state in the photon direction (O = ir) is more then 10 times less than the 
production cross section of vector (co)-state. These calculation results become similar to each other 
at high energy and in the limit as m q — > (see Fig. 2). It is seen from these figures that the (cri- 
state production in the backward direction (in the photon direction) is suppressed in comparison 
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with the production in the forward direction (the factor of m 2 (s)). As a rule at fixed target facilities 
the research of production is only possible in the photon direction. Therefore the contribution of 
perturbative recombination must be small. In addition to that the charge asymmetry should 
decrease with the increasing of beam energy. On the other hand there is no suppression at the 
production yield alongside to hadron direction (see eq. (6)). It leads to energy-independent charge 
asymmetry. 

Before proceeding further, it is necessary to remove perplexity which could appear after com- 
paring the predictions (6) and (7) with the results following from Fig. 2. Indeed one can see from 
Fig. 2 that the yields of vector and pseudoscalar mesons are practically identical in the back- 
ward production (O ~ 7r), while the vector mesons dominate in forward production (0 ~ 0). 
However, these results contradict (6) and (7). The point is that the mathematical limits of the 
cross section values evaluated in [9] are correct. But these quantities have no physical mean- 
ing. As it is shown in our study the production cross sections for vector and pseudoscalar states 
are practically the same in the kinematical region —0.99 < cos0 < —0.8, whereas the limit 
regime establishes only at cos < —0.99. The integral cross section value in the kinematical 
region — 1 < cos© < —0.99 is negligible in comparison with the integral cross section in the 
kinematic region —0.99 < cos© < —0.8. It was revealed also that in the forward production 
at 0.9 < cos© < 0.999 the vector mesons dominate: while at cos© = 1 the ratio between 
the vector meson production and the pseudoscalar production equals to 1, but at cos© = 0.999 
that ratio is about 7. The contributions from region 0.999 < cos© < 1, as well as from region 
— 1 < cos© < —0.99 into the total cross section are negligible. Now it is clear that the predictions 
(6) and (7) are beyond both the model accuracy and the experimental potentialities. Therefore 
the physically reasonable predictions at large s are as follows: 



da(7 + q (cg)( 3 Si) + c) 



da(>y + (cg)( 1 s ) + c) 
dcr(7 + q -> (cg)( 3 si) + c) 



> 1, da(>y + q-> (cq) + c)\ e 2 ~ - (8) 



2 



s 



~ 1, dcr( 7 + q -> (cq) + c)| 9siW)8>m 2 ~ ^. (9) 



da(j + q -> (cg)( 1 s ) + c 
Let us consider in detail the both perturbative recombination and photon-gluon fusion. We 
restrict ourself to the calculation of production asymmetry for D*-mesons, because we expect the 
suppression of the D-meson production in the photon direction in comparison with the _D*-meson 
production. 

For the numerical analysis we choose the kinematical cuts of COMPASS experiment, t.e.: the 
averaged energy of the muon beam is 160 GeV, the photon virtuality Q 2 < 1 GeV 2 ; xp > 0.2, 
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where xf is the Feynman variable for _D*-meson; z 7 > 0.2, where z 1 is the fraction of photon 
momentum carried out by _D*-meson. 

Fig. 3 presents the cross sections for D* + and D*~-meson production calculated in the frame 
work of the perturbative recombination mechanism has been shown as well as the contribution due 
to photon-gluon fusion. The perturbative recombination contributes to the charge asymmetric part 
of the total cross section, while the photon-gluon fusion contributes to the charge symmetric one. 
The predictions of perturbative recombination and photon-gluon fusion for D*° and Z)*°-meson 
production are presented in Fig. 4. 

In the kinematical region under consideration the perturbative recombination dominates due 
to a valence quark of the initial proton, as it is seen in Fig. 3 and 4. On the face of it, it seems 
that cr(D-) < a(D + ) and a(D°) < a(D ). However, it is worth to note that D-(D°) -meson in 
the perturbative recombination is produced in association with the c-quark (see Fig. 5 and 6). To 
estimate the charge asymmetry value resulted from from a pertubative recombination one needs 
to keep in mind that c-quark also produces the charm particles. Let us suppose, that c-quark 
hadronizes into final states with probabilities Wd* (c — > D*), Wd (c — > D), Wd s (c — > D s ) and 
W\ c (c — > A c ) (each do not necessarily equal to the appropriate probabilities at large c-quark 
momenta). To calculate the contribution of c-quark factorization process into the cross section 
one needs to take into account only particles, which are produced by c-quark in the kinematical 
region of COMPASS experiment. 

To explain the assertion stated above, let us consider an example. To estimate the yield of the 
/}*~-meson produced in the recombination process, — > D*~(cd) +c, we apply the following cuts 
xp{D*~) > 0.2, z^(D*~) > 0.2 and do not care about c-quark. To estimate the contribution of 
the remained c-quark, then we apply the cuts on xp > 0.2 and z 7 > 0.2 for the product of c-quark 
hardronization and do not care about _D*~(c<i)-meson. 

Since an interaction energy is closed to threshould one, it is reasonable to suppose that the total 
momentum of c-quark carried out by the charm meson. However, one can not exclude that the 
charm meson carries out only a part of the c-quark momentum, as it takes place at high interaction 
energies. That is why we consider two extreme cases: 

1. the total momentum of c-quark is passed to meson; 

2. An energy loss of c-quark is described by Kartvelishvili fragmentation function [11]: 

f{z)^z 2 ^{l-z D ) (10) 
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The charge asymmetry of D* -meson production is described by the equation: 



rr D * + — rr D *~ 
AD*+ _ "total "total 

~ a D*+ 7 D- • 
"total ~ "total 



Let us o"^ be the cross section value for meson production in the photon-gluon fusion; a® q is 
the cross section value of meson production in the perturbative recombination with light quark 
q; while cr° is the cross section value of the "nonrecombinated" c-quark production. Using these 
notations the asymmetry has the form as follows: 

KT + <r ) + + °* + ^ + ^) + 2 <; + ' 

For the case (1) when the meson carries out the whole c-quark momentum our calculation leads 
to the following result: 

A? * + ~ -0.03. 

For the case (2) when an energy loss of c-quark is described by Kartvelishvili fragmentation 
function [11], one gets: 

Af + ~ -0.17. 

For the neutral meson production one has another predictions: 

A? *° ~ -0.08, 
Af ° ~ -0.21. 

One can see that the asymmetry values both for neutral and charged mesons strongly depend on 
assumptions about the fragmentation mechanism for the "nonrecombinated" c-quark. Therefore a 
lack of knowledge about this mechanism does not allow to calculate asymmetry values A D * + and 
A D *° with reasonable accuracy. 

Nevertheless, the more accurate prediction really exits in the considered kinematical region. In 
spite of the absence of s-quark inside an initial nucleon the perturbative mechanism leads to the 
charged asymmetry in -D<,-meson production. As it mention above, the production of the mesons 
with c-quark in the perturbative recombination occurs with an additional production of c-quark. 
This c-quark can hadronize into Df. Therefore an additional yield of D*~ and Z)*°-mesons leads 
to an additional yield of _D+-mesons. 

The asymmetry for .D^-meson production is described by following expression: 

A Dt = -^.(°;-^ + 4- g ;) (13) 
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Our estimates are given below: 

Af s+ ~ 0.57 

and 

~ 0.64 

Thus, the perturbative recombination mechanism leads to the essential charge asymmetry in 
the process of -D s -meson production. The value of such asymmetry depends weakly on assumptions 
about the fragmentation mechanism for the "nonrecombinated" c-quark. 

3. CONCLUSION 

The perturbative recombination mechanism, which violates a factorization theorem in charmed 
hadron production causes the charge asymmetry in production of these particles. To calculate 
correctly the asymmetry value, ones need to take into account the hardonization of "nonrecombi- 
nated" charm (anticharm) quark, which is produced together with the anticharm (charm) quark 
to be recombinated with the light constituent quark of initial proton. 

In the considered process of the charm photoproduction into the region of photon fragmentation 
the nonzero value of charge asymmetry takes place due to the backward scattering of valence quark 
of the initial hadron. In the discussed kinematical region the asymmetry becomes negligible with 
energy increasing. In the condition of COMPAS experiment the asymmetry value is small in 
Z}*-meson production, but it is essential in D s -meson production. 

Therefore one can observe an interesting phenomenon: the production of meson, which 

contains the light valence quark from the initial proton, leads to the charge asymmetry in the 
production of -D s -mesons, which does not contain the valence quark from the initial proton. In 
spite of small D s meson yield, the charge asymmetry A D * is large and it depends insignificantly 
on details of a fragmentation model. Note that the first time a such effect was considered in [12] 
where a possible asymmetry in f^-meson production in pp-collisions was investigated. This process 
had been also investigated in paper [9]. 
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da /(fcos0, nb 




j i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i 
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cos 6 

Figure 1. The cross section distributions over the scattering angle cosine of (dc)-meson in 
c.m.s. for the process dj — > (dc) + c both for the vector (solid curves) and pseudoscalar (dashed 
curves) meson states. Ours predictions (bold curves) have been performed in comparison with 
the predictions of work [9] (thin curves). y/ s d~ ( = 5 GeV, m c — 1.5 GeV, rrid — 0.3 GeV. 




Figure 2. The same as in Fig. 1 but for ^fsd^ — 100 GeV, m c = 1.5 GeV, and irid = 0.1 GeV. 
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da/dp c ± , nb/GeV da/dxp, nb 
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Figure 3. The cross sections as a function of (a), x F (b), and z 1 (c) for the following subpro- 
cesses of charged charmed meson photoproduction: £>*^-meson production in the perturbative recombi- 
nation process drf — ► Z?*~ + c (solid curve); Z)* + -meson production in the perturbative recombination 
process drf — > Z?* + + c (dashed curve); the production of D* + and £>*~-mesons in photon-gluon fusion. 



13 



da/dp c ± , nb/GeV da/dxp, nb 




Figure 4. The cross sections as a function of p± (a), xf (b), and z 1 (c) for the following subpro- 
cesses of neutral charmed meson photoproduction: Z)*°-meson production in the perturbative recombi- 
nation process uj — > D*° + c (solid curve); Z?*°-meson production in the perturbative recombination 
process uj — > D*° + c (dashed curve); the production of D*° and Z)*°-mesons in photon-gluon fusion. 
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da/dp c ± , nb/GeV da/dxp, nb 




0.2 0.4 0.6 0.8 1 



Figure 5. The cross sections as a function of the kinematic variables of "nonrecombi- 

nated" charm quark pj_ (a), x c F (b), and (c) for the following subprocesses of charged 
charmed meson photoproduction: Z?*~-meson production in the perturbative recombination process 
dq — ► + c (solid curve); £>* + -meson production in the perturbative recombination process 

dj — > D*+ + c(dashed curve); the production of D* + and D*~-mesons in photon-gluon fusion. 
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da/dp c ± , nb/GeV da/dxp, nb 
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Figure 6. The cross sections as a function of the kinematic variables of "nonrecombi- 

nated" charm quark p c ± (a), x c F (b), and (c) for the following subprocesses of neutral 

charmed meson photoproduction: £)*°-meson production in the perturbative recombination process 
u-y — > D*° + c (solid curve); D*°-meson production in the perturbative recombination process 
uj — > D*° + c (dashed curve); the production of D*° and i)* -mesons in photon-gluon fusion. 



